PROCEEDINGS 

°F SCIENCE 



00 

o 
o 

(N 

> 

O 

O 



> 

cn 
(N 

d 

oo 

o 



Charm quark system in 2 + 1 flavor lattice QCD using 
the PACS-CS configurations 



PACS-CS Collaboration : Y. Namekawa*'^ S. Aoki^ S N. Ishll ', K.-l. Ishlkawa'^, 
N. Ishlzuka" ^ T. IzubuchI % D. Kadoh", K. Kanaya' , Y. KuramashI" ' , M. Okawa^, 
Y. TaniguchI" ^ A. Ukawa' ^ N. Ukita ', T. Yoshle'' '^ 

" Center for Computational Sciences, University of Tsukuba, Tsukuba, Ibaraki 305-8577, Japan 
^Graduate School of Pure and Applied Sciences, University of Tsukuba, Tsukuba, Ibaraki 
305-8571, Japan 

'^Riken BNL Research Center, Brook-haven National Laboratory, Upton, New York 11973, USA 
"^Graduate School of Sciences, Hiroshima University, Higashi-Hiroshima, Hiroshima 739-8526, 
Japan 

"Institute for Theoretical Physics, Kanazawa University, Kanazawa, Ishikawa 920-1192, Japan 

We study heavy-heavy and heavy-light quark systems for charm with a relativistic heavy quark 
action in 2 + 1 flavor lattice QCD. Configurations are generated by the PACS-CS Collaboration at 
the lattice spacing is a = 0.09 fm with the lattice size of 32^ x 64 employing the 0(a)-improved 
Wilson quark action and the Iwasaki gauge action. We present preliminary results for the char- 
monium spectrum and the D and Ds meson decay constants evaluated at 3.5 MeV< mud < 12 
MeV with TOs around the physical value. We investigate the dynamical quark mass dependences 
of the hyperfine and the orbital splittings. The decay constants are compared with the recent 
experimental values. 



The XXVI International Symposium on Lattice Field Theory 
July 14-19 2008 
Williamsburg, Virginia, USA 



* Speaker. 

^ E-mail: namekawaOccs.tsukuba.ac.jp 



© Copyright owned by the author(s) under the terms of the Creative Commons Attribution-NonCommercial-ShareAlike Licence. 



http://pos.sissa.it/ 



Charm quark system in 2 + I flavor lattice QCD using the PACS-CS conflgurations 



Y. Namekawa 



1. Introduction 

Precise determination of physical quantities for heavy quark systems provides us with an op- 
portunity to search for new physics beyond the standard model. For this purpose lattice QCD 
should be a powerful tool. However, the use of conventional lattice quark actions aie problematic 
because of large cutoff errors due to the heavy quark masses. So far several approaches to avoid 
this problem have been employed for the study of heavy quark physics on the lattice[|I]]. Our choice 
is to employ the relativistic heavy quark action of Ref. [Q]. This formalism allows us to take the 
continuum limit in which mga corrections are controlled by a smooth function. In fact the cutoff 
errors are reduced from 0{{mQa)") to 0{f{mQa){aAQCDy) where /{niQa) is an analytic function 
around niQa = 0. 



2. Simulation parameters 



We simulate the charm quark system with the relativistic heavy quark action of Ref. ^ on 
the 2+1 flavor lattice QCD configurations which are generated by the PACS-CS Collaboration 
employing the nonperturbatively 0(a) -improved Wilson quark action with = 1.715[|3|] and the 
Iwasaki gauge action. The lattice size is 32^ x 64 whose spatial extent is L = 2.9 fm with the lattice 
spacing of a = 0.09 fm. The dynamical up-down quark mass ranges from 67 MeV down to 3.5 
MeV which is close to the the physical value. The details of the configuration production and light 
hadron physics that emerges from them are described in Refs. [Q, ^, ^, ^. Table [l] summarizes the 
simulation parameters and the statistics of the configuration sets we have used for the heavy quark 
measurements. The number of the source points is to be quadrupled. We emphasize that the data 
point at (KTud, K's)=(0.137785,0.13660) is almost on the physical point: the up-down quark mass is 
only 30% heavier and the strange quark mass is almost exactly at the physical value. 

The relativistic heavy quark action proposed in Ref. [0] is given by 



Sq = Y^Q,D,,yQy, 



(2.1) 



-K"h 



^ Fij (x) Oij + C£ ^ Fi4 (x) a;4 



(2.2) 



where we are allowed to choose = 1, while the other four parameters v, r^, cb, ce should be 
adjusted in the mass dependent way. We use the one-loop perturbative values for r,, cb and ce 
evaluated in Ref. [^]. For the clover coefficients cb and ce we incorporate the nonperturbative 
contributions at the massless limit adopting the procedure cb,e = {cB,E{niQa) — cb.e{0))^ + ^sw- 
At each simulation point, the parameter v is nonperturbatively determined from the dispersion 
relation for the spin-averaged IS state of the charmonium: 



Eipf = E{pY + ci^\p\ 



(2.3) 
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Table 1: Simulation parameters. Quark masses are perturbatively renormalized in the MS scheme. The 
renormaUzation scale is ;U = 1/a for fCud < 0.137785 and IJ.—2 GeV for the physical point. 





mM/(At) [MeV] 


m^^(^) [MeV] 


#conf 

measured total/MD time 


#source 


0.13770 0.13640 


12.3(2) 


90(1) 


400 800/2000 


1 


0.13781 0.13640 


3.5(2) 


87(1) 


100 198/990 


1 


0.137785 0.13660 


3.5(1) 


73(1) 


90 200/1000 


1 


physical point 


2.53(5) 


72.7(8) 








1 1.02 1.04 1.06 1.08 1.1 1.12 1.14 1.16 1.18 1.2 
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Figure 1: Nonperturbative tuning of v with K-h = 0.11022 at (K-ud, K-s)=(0.13770,0.13640). A red symbol 
denotes the result for the nonperturbative v. 

where v is adjusted such that the effective speed of light Ceff becomes unity. Figure [l] shows an 
example of the nonperturbative tuning of v with = 0.11022 at (jCud, k:s)=(0. 13770,0. 13640). In 
order to search for the physical charm quark mass point, we employ two values of the hopping 
parameter of a heavy quark k^. The values of Kh are chosen to sandwich the physical charm quark 
mass. 

3. Heavy-heavy system 

Let us first investigate the charmonium spectrum. At each combination of (jfud, JcJ we deter- 
mine the physical point of the charm quark by the condition that the mass of the spin-averaged IS 
state reproduces the experimental value, M{\S) = (M^^ + 3M,/^)/4 = 3.0677(3) [GeV][@]. For 
this purpose we linearly interpolate the results for M{IS) at two values of the hopping parameter 
K"h. Figure ^ illustrates this procedure for the case of (K"ud, K"s)=(0. 13770,0. 13640). 

In Fig. |3| we plot results for the mass of the orbital excitation m),^^{lP) — mj/^(15) at the 
physical point of the charm quark mass. It is hard to detect the dynamical quark mass dependence 
within the errors in the range of 3.5MeV<mud<12MeV and 73MeV<ms<90MeV. A very short 
chiral extrapolation is made employing a linear function of the up-down and the strange quark 
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Physical point of charm quarl< ^ 



8.8 8.85 8.9, , 8.95 9 9.05 9.1 

1 / Ki, 

Figure 2: Interpolation of the spin-averaged 15 charmonium mass to the physical point as a function of K^. 
Errors are within symbols. 



masses: 

mx,,{lP)-mj/^{lS) = a + pm^,i + Ym,. (3.1) 

In Fig. m we find that the extrapolated value and its error are almost identical to the result at 
(K"ud, K"s)=(0. 137785,0. 13640). This illustrates how close our simulation points are to the physi- 
cal point. The result at the physical point is consistent with the experimental value within the 
error. 

Figure ^ shows the results for the hyperfine splitting mj/^ — mf,^,. As in the orbital excitation 
we find little dynamical quark mass dependence. We extrapolate the results to the physical point 
employing a linear function of the dynamical quark masses. The extrapolated value, which is 
essentially determined by the result at (JCud, k:s)=(0. 1377 85, 0.1 3640), shows a 10% deficit from the 
experimental value. Possible sources of the discrepancy are 0{g^a) effects in the relativistic heavy 
quark action, dynamical charm quark effects and disconnected loop contributions. A recent 2-1-1 
flavor lattice QCD calculation with the highly improved staggered quarks[|l^] shows a similar value 
to ours. 

We compare our results for the hyperfine splitting in A/j = 2 + 1 QCD with the previous 
Nf = 0,2 results|ll, |l^] in Fig. ^. We observe a clear trend that the results become closer to the 



experimental value as the number of the flavor is increased. The dynamical quarks give significant 
contributions to the hyperfine splitting. 



4. Heavy-light system 

For the heavy-light system we focus on the D and Dg mesons and their decay constants mea- 
sured at (K"ud, K"s)=(0. 137785,0. 13660). This data point is so close to the physical point that wiud 
corrections in the results could be smaller than the statistical errors. We employ perturbative values 
for the renormalization factor and the improvement coefficients of the axial vector current evalu- 
ated in Ref. [p3|]. For c^^ we incorporate the nonperturbative contribution at the massless limit by 

= {c^^{mQa)-c^^{0))^ + cY with = -0.03876106(0]. Figure | compares our results 
for the D and D., meson masses with the experimental values [p3|]. They are consistent within 
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Figure 3: Orbital excitation m^^,^ — nij/^ as a function of ni^^' . Vertical dotted line denotes the physical 
point. 
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Figure 4: Hyperfine splitting of the charmonium as a function of m^^'. A vertical dotted line denotes the 
physical point. 



the eiTors. It is noteworthy that the physical charm quark mass determined from the heavy-heavy 
system successfully reproduces the heavy-light meson masses. The results for the decay constants 
are shown in Fig. ^ where we also plot the recent 2-1-1 flavor lattice QCD results with relativistic 
heavy quark actions [[l^, 17] for comparison. Although we find a sizable discrepancy between our 



result and the experimental value for fjj^ , we should analyze the full configuration set and improve 
statistics before we derive any conclusions. 

In Fig. ^ we plot the ratios of fo^ to //> and fo^ to /k in which uncertainties coming from the 
perturbative renormalization factors and the lattice cutoff should cancel out. For both cases our 
results show larger values than the experimental ones, which originate from the discrepancy found 
in Fig. 0. 
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Figure 5: Comparison of hyperfine splittings in Nf = 0[pl]], 2jl^ and 2 + 1, together with the experimental 
value. All the lattice results are obtained at a^' « 2 GeV. 



1.90 
1.88 
1.86 
1.84 
— 1.82 
^ 1.80 

1.76 
1.74 
1.72 

1.70 



J3=1.90, 32^x64 




. >K 4 










Experiment ^ 






HPQCD and UKQCD-^ 






k;v=0. 13660^ 



> 

C5 1-97 



'/3=1.90, 32^64 














Experiment ^ 






HPQCD and UKQCD ^ 






ic,=0. 13660 



Figure 6: D (left) and Ds (right) meson masses as a function of ni^]^^- A vertical dotted line denotes the 
physical point. 
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